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1. Introduction 

The Inert Doublet Model (IDM) was introduced to accommodate or explain neutrino 
masses [1] and independently, to alleviate the little hierarchy problem while also providing 
a dark matter (DM) candidate [2] . The model represents a very minimal extension of the 
Standard Model (SM), it just contains an extra weak scalar doublet, which is odd under 
an unbroken Z2 symmetry, rendering the lightest member stable. The other members of 
this doublet are another neutral particle and a pair of charged ones. These particles can 
all be produced at colliders via their couplings to clectroweak gauge bosons, subject to the 
constraint of the Zi symmetry. The collider phenomenology has been explored in [3,4] and 
the Early Universe phenomenology has been studied in some detail in [5] and [6] . 

While the IDM has many attractive features, simplicity being an important one, it was 
felt that the introduction of CP violation in the scalar sector would make the model more 
attractive, therefore an extension to a Two-Higgs-Doublet Model (2HDM) plus an inert 
doublet model was proposed [7]. This also allows for an alleviation of the little hierarchy 
problem. We shall refer to the resulting model as IDM2. 

It has been found that the IDM permits a DM particle with a mass in one of three 
regions: light (m ^ mw) [11,12], medium (m ~ my^') [2,5] or heavy [m <; 535 GeV) [5,6,8]. 
Two of these mass regions (m ^ mw and m ^ 500 GeV) were also found to yield solutions 
for the IDM2 [7]. The aim of the present paper is to explore the IDM2 in more detail, 
determine the allowed mass regions for the DM particle, its dominant annihilation channels, 
and the corresponding mass regions for the lightest Higgs boson Hi. Furthermore, we 
will confront the model with constraints from direct-detection experiments, and briefly 
comment on possible signals in LHC experiments. 

The paper is organized as follows. In Sec. 2 we review the model, and in Sec. 3 
we discuss the theoretical and experimental constraints. In Sec. 4 we consider various 
annihilation channels and in Sec. 5 we present the scan strategy adopted to search for 
allowed regions in the parameter space. In Sec. 6 we give an overview of allowed regions 
of DM particle masses, whereas in Sees. 7 and 8 we explore in more detail parameters that 
are compatible with all the constraints in the low and high DM mass regions, respectively. 
Then, in Sec. 10 we discuss constraints from direct detection experiments, in Sec. 11 we 
briefly explore relevant LHC phenomenology, and in Sec. 12 we summarize. The relevant 
couplings of the model are given in an appendix. 

2. IDM2 — model and notation 

2.1 Fields and potential 

We denote the doublets of the 2HDM as 




(2.1) 
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where = vf + and tan^ = V2/V1. The inert doublet is decomposed as 



V 



iA)/V2 J 



(2.2) 



it transforms under an unbroken Z2 symmetry as 77 ^ —r] which ensures that tj couples 
only bilinear ly to other scalars and to the gauge sector. All other fields remain neutral 
under this transformation. 
The potential reads 



y($l,$2,??) = 142($1,$2) + ^3(r?) + Vi23($l,$2,r?) 

where the 2HDM and inert-sector potentials read 

^12(^1,^2) = |m?i$|$i +m|2$2^2 + 



(2.3) 



+ 



1$2 + h.c] } 

^($t$^)2 + ^($t$2)' + Xs{^l^lM^2) 
+ A4($l$2)($|$l) + ^ [A5($I$2)' + h.c] , 

whereas their mutual couplings, bilinear in the Z2-odd field r/, are given by 

^123(^1, $2, ??) = Aii33($l$i)(r/1"r?) + A2233(^2^2)(??^^/) 
+ Ai33l(^l??)(??"^^l) + A2332(^k)(^'^^2) 

' 'A2323(^k)^ + h-C. 



(2.4) 
(2.5) 



Al313(^ir?)^ + h.C 



(2.6) 



Here, A1133, A2233, A1331 and A2332 are real, whereas A1313 and A2323 can be complex. It has 
been assumed that the total Lagrangian is symmetric under Z2 which acts as $1 — > — $1 
and ur —7- —un (all other fields are neutral) and therefore the FCNC are eliminated at 
the tree level. In order to allow for CP-violation in the scalar sector we allow for soft Z'2 
breaking by the mass 777,^2^1^2 + h.c. (so that the renormalizability is preserved). Note 
that, as a consequence of the unbroken Z2 associated with the inert doublet, there is no 
mixing in mass terms between $1^2 and r] [7]. 

The model we discuss here bears some similarity to the three-Higgs-doublet Weinberg 
model [9,10]. The Weinberg model was constructed as the minimal model that accom- 
modates CP violation in the scalar sector together with natural flavour conservation. In 
fact, the symmetries of the Weinberg model are exactly the same as those we impose here: 
Z2 X Z2, where Z2 is responsible for flavour conserving Yukawa couplings of $1 and $2 
while Z2 decouples 77 from fermions. One difference between the Weinberg model and the 
present one is that we allow for soft Z2 violation by the mass term 7nf2^i^2 + h.c. The 
second is that, in the Weinberg model, all three doublets develop non-zero vacuum expec- 
tation values, while in our case {rj) = 0. That implies that there is no mixing in the mass 
matrix between $1^2 and 77; in particular, 77^ decouples from and H^, so that there is 
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no CP-violation mediated by charged scalars. However, since 7^ 0) CP is violated in 
the neutral non-inert scalar sector in the same way as in the 2HDM. 

In the absence of the potential (2.6), one would have to require > in order to 
ensure (rj) = 0. However, the non-zero expectation values of the other fields, vi/\/2 and 
V2/V^, lead to an overall coefficient of the term that is bilinear in 77 of the form 





( Vi 




= r\v2 







fn^ + (A1133 + A1331 -I- Re Ai3i3)-^ -I- (A2233 + A2332 + ReA2323)^- (2-7) 
This must be positive, but is not necessarily positive. 
2.2 Mass eigenstates of the IDM2 

The neutral states of $i,$2 will in general mix to form three neutral states Hi, H2, H^. 
These are linear combinations of 771, r]2, and 773, 



(2.8) 



where 773 = — sin^xi-|-cos Px2 is orthogonal to the neutral Goldstone boson G° = cos Pxi + 
sin f3x2 and the rotation matrix R is parametrized in terms of three angles ai , a2 and 
according to the convention of [13]. 

For the quartic couplings describing the interaction between 77 and $1 and $2, we 
adopt for simplicity the "dark democracy" : 

Aa = A1133 = A2233) 
Aft = A1331 = A2332, 

Ac = A1313 = A2323 (real), (2.9) 

then V123 is invariant under $1 o $2- The dark-sector masses can be written as: 

=m2 + iAa^;^ 
M^ = ml + i(Aa + Ab + = Mj± + i(Afe + X^y, 
Ml = ml + i(A„ + Xh- Xcy = + i(A6 - X^y, (2.10) 

where ttt^ is a mass parameter of the 77 potential (2.5). We shall take the scalar, S, to be 
the DM particle, i.e., Ms < Ma- The other choice would simply correspond to Ac — >■ —Ac, 
without any modification of the phenomenology described in the following. 
It is instructive to invert the relations (2.10): 

Aa = ^ (m2± - rri^) , (2.11a) 
Ab = ^ (m| + Ml - 2M2±) , (2.11b) 
^c = 4(M|-Ml). (2.11c) 
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Thus, these couplings of the inert doublet to the non-inert Higgs sector can be expressed 
in terms of the mass splittings (including also the soft mass parameter m,,). 
It is convenient to introduce the abbreviation 

\L=\{\a + h + K)= \^ \ (2.12) 

From Eq. (2.7), the condition (r/) = can now be written as 

ml + Xlv'^ = MI> 0, (2.13) 
which is automatically satisfied by our choice of input parameters. 



3. Theoretical and experimental constraints 

We here present a summary of the constraints imposed on the model. Some of the the- 
oretical ones (positivity, in particular) are absolute^ whereas the experimental ones are 
quantitative, and subject to experimental precision. 

3.1 Theoretical constraints 

• CP violation 

We do not impose CP conservation on the neutral Higgs sector. The amount of CP 
violation that remains after all constraints are imposed is determined afterwards. 
For a detailed discussion of the conditions for CP to be violated in this model, see 
Appendix B of ref. [7]. 

• Stability or positivity 

The potential should be bounded from below for any values of the fields $i, $2 and 
r/. This condition is rather involved for the potential (2.3). The full set of conditions 
are given in Appendix A of [7]. For the somewhat simpler case of dark democracy 
considered here, we must impose 

Ai > 0, A2 > 0, A^ > 0, (3.1) 
Ax > — \/AiA2, \y > \y > (3.2) 

Aj, > V (a^Ax - A^ > -^(A^Ai - A2)(A^A2 - A2)) , (3.3) 

where 

Ax = A3 + min(0,A4-|A5|), (3.4a) 

Aj; = Aa + min(0,Ab-|Ac|). (3.4b) 

With Ms < min(ikOi,M^±), we have 

\ = ^{^s-ml). (3.5) 



-5- 



In [7] we constrained the potential further, by requiring V12, V3 and V123 individually 
to satisfy positivity. The condition (3.3) was then replaced by 



Aa > max(0, -2A6, -A^ ± Ac). 



(3.6) 



In terms of masses, this means for the two cases: 



Ms<Ma< M^± : ml + - < M|, 
Ms < M^± <Ma: ml< M|. 



(3.7a) 
(3.7b) 



In the present study we go beyond the domain of parameters allowed by (3.6) or, in 
terms of masses, (3.7). The latter condition allows for checking positivity "once and 
for all" , for a given set of inert-sector parameters. The full condition (3.3) depends 
also on the parameters of the non-inert sector, and must thus be checked for each 
point. The "reward" is that less parameter space will be excluded. 

• Electroweak symmetry breaking 

In order to break the electroweak symmetry spontaneously, the vacuum expectation 
values of $1 and $2 should be non-zero. The most general form of the vacuum can 
always be written in the form 



where f 1 > and u, V2 and 6 are real numbers. Non-zero u would imply spontaneous 
^^(1)em violation, so electric charge non-conservation. However, it has been shown 
in [15] that if a local charge-conserving minimum exists, then there can be no charge- 
breaking minima (there may exist a stationary point with u 7^ 0, which is a saddle 
point). Therefore from here on we assume u = 0. Nevertheless, the potential of 
the 2HDM can have more than one charge conserving minimum. In that case it is 
important to make sure that the theory is expanded around the global one, so that 
the issue of tunneling to the lowest one does not appear. In our approach we start 
out by assuming 6 = 0^ and choosing a value of tan/3. Then the parameters of the 
potential are adjusted so that vi = vcosP and V2 = f sin/3 satisfy the stationary- 
point equations. Since all scalar mass squares that we consider are positive the 
stationary point must be a minimum. However wc do not beforehand know if our 
starting minimum is the global minimum of the potential since the 2HDM allows for 
more than one minimum. In our scans, we look for points ("good") in the parameters 
space that satisfy all theoretical and experimental constraints. For each "good" point 
we also check if our starting minimum (vi, ^2) is the global minimum. If there exists 
a deeper minimum we discard the point. Of all the "good" points that satisfied all 
other restrictions, approximately 7% were thus discarded because a deeper global 

^Note that for a given minima it is always possible to make its location real by performing an appropriate 
global phase rotation. 




(3.8) 
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minimum exists. Note that if the Universe was indeed in a state corresponding to a 
false vacuum, then the tunneling to the true vacuum would in principle be possible. 
If the tunneling time was shorter than the Universe age that could have important 
cosmological consequences. This, however lies beyond the scope of the present study. 
Therefore we restrict ourselves to the case of global minimum only. Note that the 
presence of the third doublet 77 does not influence the above arguments since (r/) = 0. 

• Unitarity and perturbativity 

We impose unitarity on the non-inert Higgs-Higgs-scattering sector [16-18]. At large 
values of tan the soft mass parameter /x is rather constrained, /x ~ M2 ~ Mfj± , as 
discussed in Refs. [19,20]. Furthermore, perturbativity is imposed, in the form 

Xi, |aj|, |aj|, ^^cot/5, Aa,A6, Ac < 47r. (3.9) 

V V ■\/2v 

Here, aj and aj are coefficients of the CP-even and odd parts of the Yukawa couplings 
[21]. For the couplings to t {b) quarks, they are given by aj = Rj2/ sin/3 and aj = 
—Rjz/ tan^ {cLj = Rji/ cos (3 and aj = —Rj^tanP). 

• The little hierarchy 

In order to avoid excessive computational requirements, and in distinction from the 
approach of [7], we will not a priori impose an alleviation of the little hierarchy. 
However, it turns out that for Ms 2± 75 GeV, the masses of the non-inert Higgs 
sector can be lifted to rather high values, as will be discussed in Sec. 7. This provides 
for a considerable alleviation of the little hierarchy problem around Ms 75 GeV. 

3.2 Experimental constraints 

We impose a variety of relevant experimental constraints. These can be grouped as follows: 

• Charged-Higgs sector 

The charged-Higgs sector is constrained by several observablcs. The B — > Xg^ data 
constrain low values of tan /3 and low Mjj± , the details of which depend on QCD 
effects [22-25]. Likewise, the Bq — Bq mixing constrains low values of tan/3 and low 
Mff± [26-28], whereas B Drvr and B tUt-X constrain low values of Mjj± and 
high values of tan /3 [29-34]. 

• Neutral-Higgs sector 

The LEP2 Higgs boson searches have given limits on the coupling of the lightest 
Higgs to the Z and to bh [35]. At low values of tan^, and low Mff±, the well- 
measured T{Z — > bb) decay rate also constrains charged-Higgs contributions, and 
to a much lesser extent, neutral-Higgs couplings [36,37]. There arc also bounds 
stemming from the clcctroweak precision data, in terms of the T and S parameters 
[38-40]. Among these, the most serious one is the constraint on T, which basically 
is a constraint on the mass splitting of pairs of scalars: A pair of neutral scalars 
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or a pair of charged scalars of different masses botli give a positive contribution 
AT{Mf , Mj) > 0, whereas a neutral-charged pair gives a negative contribution, 
AT{M^,Mj) < 0. Such contributions must roughly cancel, in order not to violate 
the electroweak precision data. 

Furthermore, at large values of tan P, the model is constrained by the electron electric 
dipole moment [41-43], for which we adopt the bound: 

|de| < 1 X 10-27[ecm], (3.10) 

at the l-fj level. This is calculated directly from the neutral-Higgs-sector mixing 
matrix [7,42]. The muon anomalous magnetic moment [37,43,44] has however little 
impact, since the large tan P region tends to be excluded by the unitarity constraint 
[37]. For the relevant loop calculations, we use the LoopTools package [45,46]. 

• Inert-sector constraints 

The amount of dark matter has now been measured to an impressive precision [47] 

J^dm/i^ = 0.1131 ± 0.0034 (3.11) 

We estimate the model prediction of the amount of dark matter from an implemen- 
tation of micrOMEGAs [48,49]. 

For the heavier, neutral member of the inert sector, we adopt the bound obtained 
from a re-analysis of LEP data [4], approximated as Ma > 110 GeV. For the charged 
member, we adopt the LEP bound on the chargino mass [50], M^± > 70 GeV. 
This is slightly more conservative than the bound on charged Higgs bosons, M^± > 
79.3 GeV, adopted by Ref. [5]. 

Since some of these constraints are correlated, we do not accumulate their "penalties" 
in the form of an overall measure, but rather demand that each of them be satisfied to 
within 2a. 

4. Annihilation mechanisms 

In order not to over-produce dark matter in the Early Universe, annihilation channels must 
be kinematically open. These are of different kinds, depending on the mass scales involved. 
"External" reference mass scales are the W and the lightest Higgs mass scales, mw and 
Ml. If the DM mass is low compared to mw it will annihilate via the lightest (but off- 
shell) Higgs particle which then decays to bb or cc. If the mass is comparable to mw, it 
can pair-annihilate to W^W~ or ZZ. In the higher mass range, the neutral and charged 
members of the inert doublet will be ncar-dcgcncrate, and several channels will be open. 
We here review the different annihilation mechanisms, relevant in different mass ranges. 
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4.1 DM couplings 

The gauge and scalar couplings involving inert-sector fields are collected in Appendix A. 
Here we list some of the most relevant ones. 

The DM particles can annihilate via the gauge coupling: 

SSW+W- : (4.1a) 

SSZZ : , (4.1b) 

2 cos^ &w 

or to non-inert scalars via the following trilincar or quartic couplings: 

SSHj : -2iFssj^LV, with Fssj = cos l3Rji + suil5Rj2, (4.2a) 

SSHjHj : - 2i{\L - XcR%), (4.2b) 

SSHjHk : 2i\cRj3Rk3, (4.2c) 

SSH+H- : - iXa, (4.2d) 

where is defined in Eq. (2.12) and the pre-factor in (4.2a) satisfies \Fssj\ < 1, since 
R is unitary. In particular, Fssi = cos(/3 — ai)cosa2- Also, we note that the splitting 
MJ — controls the strength of the important trilinear coupling to a neutral Higgs 
field. Likewise, the couplings A^ and Ac are related to mass splittings, Xa ~ i^^^± ~ ''^^ri) 
and Ac ~ (M| - M^), see Eq. (2.11). 



4.2 Representative branching ratios 

In the medium-mass region, the early-universe abundance is controlled by SS annihilation, 
typically to bb or W~^W~. In the high-mass region, in part due to the high degree of mass 
degeneracy of the different states, there are also significant losses due to annihilations of A 
and r]^. The different losses are provided by micrOMEGAs as fractions of l/{^r>Mh^) (see 
ref . [48] ) . We refer to this variable as "loss channel" . Some representative values are given 
below, where we consider separately the low- and medium-mass region, and the high-mass 
region. 

4.2.1 Low- and medium- mass region 

In the low and medium-mass region, the annihilations mostly proceed via an intermediate 

Higgs boson, or via a W~^W~ pair or a ZZ pair: 

• For Ms = 40 GeV and Mi = 120 GeV, representative loss channels are: 

SS^l'' (4.3) 
[cc (1-9%), 

• For Ms = 80 GeV and Mi = 120 GeV, representative loss channels are: 

,,_^IW-W- (19-96%), 
[bb (2-91%), 

where the ranges relate to the scans over M^±, m^, tan/3, M^± and the a^, defined in 
Sec. 5. 
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4.2.2 High-mass region 

The inert scalar masses are generated by m^, vi and V2 such that contributes universally 
to all the masses while the splitting between them is controlled by the interaction terms 
(Aa, Xb and Ac) between the 2HDM and the inert sector. As the DM mass Ms increases, 
annihilation channels to pairs of vector bosons opens (with fixed gauge coupling constant 
strength). Therefore annihilation through intermediate Hi must be suppressed in order not 
to introduce a too small value of flnu, that implies small Aa, A5 and Ac. In other words, 
the inert sector masses must be similar. For some representative mass parameters, loss 
channels arc given below: 

• For Ms = 550 GeV and Mi = 120 GeV, representative loss channels are: 

jw+W- (14-17%), 
\ZZ (11-14%), 



ss 



w+w- 

77 

'w+w- 

{ZZ 



AA 



(14-17%), 
(13 - 16%), 

(5 - 6%), 

(8-9%), 
(6 - 7%), 
(6- 



• 



• For Ms = 3000 GeV and Mi = 120 GeV, representative loss channels are: 

'W+W- (4-12%), 
n+n- 

SS- 



H+H- 


(3- 


- 11%), 


ZZ 


(1- 


-4%), 




(1- 


- 3%), each j. 


'w+w- 


(1 


- 9%), 


H+H- 

< 


(1 


- 9%), 


ZZ 


(2 


-7%), 


HjHj 


(2 


- 7%), each j 



AA 



w+w- 

H+H- 
ZZ 



HjHj 



Sri^,Arf 



ZW^ 

HjW^ 

HjH^ 



(3 - 10%), 

(3-9%), 

(1-2%), 

(1-3%), eachj, 
(1-9%), 

(1-10%), eachj, 
(1-10%), eachj. 



(4.5a) 



(4.5b) 

(4.5c) 
(4.5d) 



(4.6a) 



(4.6b) 



(4.6c) 



(4.6d) 
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where again the ranges relate to the scans over M^±, m^, tan/3, Mh± and the a^, defined 
in Sec. 5. 

5. Parameters and Scan Strategy 

5.1 Model parameters 

The model contains a total of 13 parameters defining the spectrum and the couplings. 
Among these, the inert-sector self-coupling. A,,, plays no role and is kept fixed. We choose 
the remaining ones to be: 

1. Ms, Ml (lowest physical masses of inert and 2HDM sectors, fixed) 

2. Mji, M^± (inert sector, physical masses, fixed). In the high Ms regime {Ms > 
500 GeV), Ma and M^± should be rather close to Ms, in order to prevent and 
Ac from becoming large (and thus lead to too much annihilation of DM in the Early 
Universe). In the low Ms regime {Ms < 100 GeV), a wider range of values is possible, 
but LEP2 data constrain Ma ^ IW GcV [4]. For M^±, we impose the constraint 
M^± > 70 GeV, adopted from the LEP searches for charginos [50]. 

3. M2, jU (2HDM sector parameters) 

4. rrirj (inert sector, soft mass parameter, fixed). Then \a,Xb,\c are all fixed. In the 
high Ms regime {Ms > 500 GeV), nirj should be rather close to Ms, in order to 
prevent A^, from becoming large (and thus lead to too much annihilation of DM in 
the Early Universe). In the low Ms regime {Ms < 100 GeV), a wider range of values 
is possible. 

5. tan^, Mh± (2HDM sector). We allow for the following variation: 0.5 < tan^ < 50 
and 300 GeV < M^± < 700 GeV. We consider a logarithmic grid in tan P, and linear 
in Mff± , typically 30 points in each parameter. Representative, allowed regions are 
shown in Sec. 9. 

6. ai, a2, 0:3 (2HDM sector). The allowed range of variation is — 7r/2 < ai^2 < 7r/2, 
and < 0:3 < 7r/2. A random set of 1000 points in this three-dimensional space is 
typically adequate. Representative, allowed regions are shown in Sec. 9. 

Prom this input, the value of M3 and all A's of the 2HDM can be reconstructed [14]. For 
the inert sector, we take A^ = 0.2 (in the notation of the IDM [2, 5] A^ = 2A2). This 
parameter has little influence on the model [5]. Prom the chosen mass input, together with 
A^, the Aa, Aft and Ac of Eq. (2.9) can be determined. 

5.2 General scanning strategy 

We scan over the parameters in a hierarchical fashion. In the notation of Sec. 5.1, we hold 
parameter fixed while scanning over pi. Por each parameter pi, we have two options: 
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(i) if we find an acceptable solution, proceed to the next value at the higher level, Pi-i, 
or 

(ii) exhaust a predefined range from to p™^. 

The latter condition is adopted if we are interested in determining the allowed range of pi. 

Scanning over this large number of parameters wc focus on regions where the 2HDM 
is known to be consistent [19]. More details are given in Sees. 7 and 8. 

It should be stressed that adopting the strategy described above, some solutions could 
be missed. However, as we do find interesting regions of parameter space that are allowed, 
we do not insist on exhausting the whole parameter space what would imply a dramatic 
increase of the CPU running time. 

5.3 Positivity and unitarity 

For the IDM, it was found [5] that certain domains in the Ms-m^ plane (denoted 
and ii2 in [5]) are forbidden by positivity. For the present model, because we have more 
parameters (in particular, the Oj parameters), we have not found any such domain excluded 
by positivity. However, if we restrict the scan to the CP-conserving limit 02 — > 0, 0:3 ^ 0, 
then some parts of parameter space are actually excluded, in particular for M^± < Ms- But 
that region is of course not interesting, since we want the DM candidate to be electrically 
neutral. 

If we impose also unitarity in addition to positivity, the allowed parameter space in 
M^±-mri starts to shrink, as indicated in Fig. 1 for Mg = 550 GeV and Mi = 120 GeV. 
For the lower range of M5- values, we do not find any such forbidden region within a 
"reasonable" range of parameters. For example, with {Ms^Ma) = (75,110) GeV, there 
is no forbidden region within 70 GeV < M,,± < 150 GeV and GeV < < 150 GeV. 
Imposing next the experimental constraints discussed in Sec. 3.2, we find rather dramatic 
reductions of the allowed parameter space, as will be discussed in Sees. 7 and 8. 

6. Overview 

Before going into a detailed discussion of allowed parameter regions, we here give a brief 
summary, comparing with the IDM, for which three mass regions were established: light 
{Ms < mw) [11,12], medium {Ms ~ mw) [2,5] or heavy {Ms ^ 535 GeV) [5,6,8]. A 
recent analysis finds a further substructure of those regions [51]. Finally, a very recent 
study [52] reports a new viable region for masses in the range mw ^ Ms ^ 150 GeV. 

We found two regions of allowed Ms values, a low-to-medium region, 5 GeV ^ Ms ^ 
100 GeV and a high region, Ms ^ 545 GeV. The allowed M^-ranges of these two models 
axe compared in Fig. 2. The other parameters axe chosen such that the allowed regions are 
maximized. In the low-to- medium DM-mass region, Early-Universe annihilation via the 
lightest neutral Higgs boson. Hi plays an important role for obtaining the observed value of 
JIdm- In the high DM-mass region, it is instead annihilation to two on-shell gauge bosons 
or two Higgs bosons that provide the correct amount of DM. In the forbidden intermediate 
range of Ms, ^dm is too low. While the present model has a few more parameters than 
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Figure 1: Regions in the Mjj± — plane that are allowed by positivity and unitarity, for DM 
mass Ms = 550 GeV and lightest Higgs mass Ah = 120 GeV. 

the IDM, we do not find solutions at as low values of Ms as were found for the IDM. A 
main restriction on the model at low Ms is the more tight up-to-date constraint on 0dm 
imposed here, the value comes out too high. In addition, our constraint on M^ is different. 
Likewise, the new viable region of the IDM [52] is only partly reproduced in our model, 
again mainly because of our more tight constraint on JIdm- 

7. Low— Medium DM Mass Regime 

For a range of DM-masses, Ms < 0(100 GeV), we have explored the range of light- 
est neutral-Higgs-boson masses, Mi > 120 GeV, for which we find consistent solutions. 
Selected results are shown in Figs. 3-6. In these figures we display, for a given set of 
(Ms, Ml) values, the allowed region(s) in the M^±-m,, plane, obtained by a scan over the 
ranges 70 GeV < M^± < 150 GeV and < rrir, < 160 GeV. 

We note that for a given set of inert-sector masses, {Ms, Ma, M^±), one may think 
of rrirj (denoted /U in [2,4,5], and not to be confused with the parameter normally denoted 
IJL in the 2HDM) as basically determining the trilinear coupling SSHi, see Eqs. (4.2a) and 
(2.12). Thus, the vertical axis in these plots is a measure of how strongly the DM particles 
annihilate via a virtual Higgs. Indeed, along the right-hand edge of these plots, we also 
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Figure 2: Top panel: Allowed Ms ranges in the IDM [2,5,6,8,11,12,52]. Bottom: Allowed Ms 
ranges in the present model, IDM2. Corresponding Mi ranges are given in Sees. 7 and 8. 

indicate some values of A^,. 

7.1 Scanning strategy 

For this low-to- medium range of Ms, we start out with a fixed value of Ma = 110 GeV 
(approximately the lower limit compatible with LEP data [4]), and then scan over M^± 
and as indicated in Sec. 5.2. As mentioned above, for M^±, we impose the constraint 
M^± > 70 GeV, adopted from the LEP searches for charginos [50]. No particular hierarchy 
is assumed, we may have Ms < Ma < or Ms < M^± < Ma- If no solution is found 
for Ma = 110 GeV, the scan is repeated for Ma = 115 GeV. If still no solution is found, 
we declare there to be no solution for the chosen set {Ms, Mi). (Only in one case was a 
solution found for Ma = 115 GeV and nothing for 110 GeV.) 

For fixed Ma, the scanning over M2 and fi is organized as follows. For Mi < 300 GeV, 
we first consider M2 = 300 GeV and /i = 200 GeV. If nothing is found, we increment /U 
to 350 GeV and 500 GeV. If still nothing is found, we increment M2 to 400 and 500 GeV, 
repeating the values of /U. For higher values of Mi, correspondingly higher values of M2 
and n are adopted, but with less splitting with respect to Mi. This choice is inspired by 
the knowledge of allowed regions in the 2HDM parameter space [19,20]. For example, for 
Ml = 300 GeV, we take M2 = 350 GeV, 425 GeV, and 500 GeV. 

7.2 Results for A/5 < 100 GeV 

Let us now comment on the main features of the plots showing allowed regions in the M^±- 
rriri plane. As pointed out above, for fixed masses of the inert sector, {Ms, Ma, M^±), the 
"soft" parameter m^, which represents the bilinear coupling in the inert sector, see Eq. (2.5), 
will also represent the trilinear couplings between the inert sector and the non-inert one, as 
expressed by Eqs. (A. 7), (2.12) and (2.11). Explicitly, the trilinear coupling SSHj vanishes 
in the limit — ?■ Ms, but can become large when either m,j <C Ms or ^ M5. 

In the case of a low- to medium-mass DM particle {Ms ^ 100 GeV), annihilations via 
a virtual Higgs boson play an important role for obtaining the correct DM density. We 
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Figure 3: Allowed regions (dark blue) in the M^± — rn^ plane, for DM mass Ms ~ 6 GeV and 
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rriri = Ms (or A^, =0), whereas the dashed line at M^± — 70 GeV (labelled "LEP bound") gives 
the adopted experimental bound. The right-most dashed line gives M,j± =110 GeV, the default 
value for Ma- 



recall that this coupling is proportional to Al, which in turn is proportional to M| — m^. 
This feature is the same as for the IDM [5], and has implications for whether or not the 
parameter region ~ Mg is allowed. However, in distinction from the IDM, here the 
coupling also contains a factor Fj, satisfying < < 1 (see Appendix A), which depends 
on tan /3 and the mixing angles of the neutral Higgs sector. Thus, for a given set of inert- 
sector parameters, one can always find non-inert-sector parameters for which this coupling 
is turned off (but positivity and unitarity constraints may make such points uninteresting). 

If the DM particle is very light, 2Ms < Mi, the intermediate Higgs particle will be 
off mass shell, and a heavier one (larger value of Mi) will be less efficient in mediating 
annihilations, because of the propagator suppression. In order to have sufficient Early- 
Universe annihilation, the S'S'lfi-coupling Al must therefore be stronger, the higher Mi 
is. Indeed, for the very lowest values of Ms, we do not find solutions for Mi = 120 GeV, 
whereas Mi = 90 GeV and a large value of give acceptable solutions (see Fig. 3). This 
is further illustrated in Fig. 4, for Ms = 20 GeV and 40 GeV, where we see that indeed 
|MJ — m^l oc |Al| must increase with Mi. 

In this region, the cut-off towards higher values of Mi is due to the AT and JIdm 
constraints. For example, at the edge of the forbidden region, for Ms = 40 GeV and 
Ml = 220 GeV, more solutions are found if we either disregard the AT constraint, or 
allow a higher value of f^DM- For this case of Mi = 220 GeV, because of the mentioned 
propagator suppression, one must have ^110 GeV. 

As Ms is further increased, two things happen. First, as Ms approaches the mass 
of the W, annihilation (in the Early Universe) via off-shell W^s and Z's starts to play an 
important role, as illustrated quantitatively by Eq. (4.4). Secondly, the lightest neutral 
Higgs can be produced resonantly via SS annihilation. This is illustrated in Fig. 5, where 
we consider Ms = 60 GeV and (left panel) Mi = 120 GeV. Only small values of the 
trilinear coupling are allowed, refiected in the plane being populated by allowed solutions 
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Figure 4: Allowed regions in the M^± — plane, for DM mass Ms — 20 GeV and Ms = 40 GeV, 
with lightest Higgs mass Mi = 120 GeV and 150 GeV. 
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Figure 5: Allowed regions in the M^± — plane, for DM mass Ms — 60 GeV, with lightest 
Higgs mass Mi = 120 GeV and 200 GeV. 



around rrirj = Ms- (The "holes" are presumably due to the "small" number of points being 
scanned over.) For a larger value of Mi, Hi is no longer produced resonantly, and a certain 
minimum value of Xl (or, equivalently, — M||) is required, as shown in the right panel 
of Fig. 5 for Ml = 200 GeV. 

As Ms reaches a value around Myy or Mz, annihilation becomes very easy via the WW 
and ZZ channels (the strength is given by the gauge coupling). However, to reproduce the 
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Figure 6: Allowed regions (dark blue) in the Af^± — plane, for DM mass Ms = 75 GeV, with 
lightest Higgs mass Mi = 120 GeV, 150 GeV, 200 GeV and 400 GeV. 



Allowed regions 
Mg: 80 GeV M,: 120 GeV 



o 



0.2 
0.1 


0.05 

0.1 



100 150 
M,. [GeV] 



Allowed regions 
Mg: 80 GeV M,: 150 GeV 



O 



0.2 
0.1 


0.05 

0.1 



100 150 
M,. [GeV] 



Allowed regions 
Mg: 80 GeV M,: 155 GeV 



-0.2 
1-0.1 


0.05 

0.1 



100 150 
M„. [GeV] 



Figure 7: Allowed regions in the M^± — m,, plane, for DM mass Ms = 80 GeV, with lightest 
Higgs mass Mi = 120, 150 GeV, and 155 GeV. 



correct DM abundance some annihilation must also take place via the Higgs (dominantly 
Hi) intermediate state. It is worth noticing that in order to obtain substantial contribution 
from that process, an increase of Mi must be compensated by larger A^^, therefore for a 
given Ml only suitable values of and non-inert-sector parameters (tan/3, Rij) provide 
an appropriate 0dm- Some such cases are shown in Fig. 6. Of course, the annihilation 
via a neutral Higgs can be resonant, in which case has to be small, as illustrated for 
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2Ms = Ml = 150 GeV in the upper right panel in Fig. 6. 

In this region of Ms ^ M]^, there is also another issue. As recently pointed out [53], the 
annihilation could be enhanced by going via WW*, where W* denotes a virtual W. Thus, 
three-body decay will set in below the two-body WW threshold. However, we do not find 
an effect as large as reported in [53]. In this region, there are two important mechanisms: 
the four-point SSWW gauge coupling, and the s-channel annihilation via an intermediate 
Higgs, the strength of which is given by Al (and thus tunable via the model parameters). 
If the intermediate Higgs is close to its mass shell, the s-channel becomes more relevant. 
However in this region of the parameter space special care is required, as the Higgs reso- 
nance contribution is very sensitive to the proper treatment of the Higgs width. This effect 
may be responsible for the reduced (as compared to [53] ) effect caused by three-body final 
states that we have found. As an example showing the importance of the three-body final 
states, wc consider one of the "good" points, {Ms,MA,Mj^±,mr^) = (75,110,86,0) GeV, 
with (Mi,M2,/i,Mj^±) = (140,300,200,389.7) GeV and tan/3 = 0.689. Allowing a virtual 
W, we find a drop of JIdm from values around 0.11 to below 0.07. If we turn off the 
s-channel, with = 0, or = Ms, the corresponding values are 0.13 and 0.05 (a larger 
effect). In view of the excessive computational requirements, we have not pursued this. 

For Ms approaching M^±, as happens near the lower bound on M^± (70 GeV), the S 
and 77^ number densities in the early universe would have been similar. In this parameter 
region, the quartic Srf'^W^j gauge coupling becomes relevant, since the Wj channel is 
kinematically open. For parameters in this region, there would be too much co-annihilation 
Srj^ —7- W^^, and JIdm would be too low. Thus, values of M^± close to Ms are not 
allowed. This is reflected as a forbidden band in Figs. 6 and 7. For higher vahics of A£,^± 
(with respect to Ms) this is not a problem, because of the Boltzmann suppression of the 
r]^ number density. 

For Ms = 75 GeV and Mi = 400 GeV, there is a strip of allowed parameters around 
M^± ~ 85 GeV, for almost all values of (see Fig. 6, lower right panel). Along this 
strip, co-annihilation of Sry^ W^j, as well as Sr)"^ — >■ W^* — >■ ud, du, cs, sc play 
an important role. On the other hand, with such a heavy Higgs boson, the role of an 
intermediate Hi is much reduced. Thus, for somewhat higher values of M,^± , the co- 
annihilation is no longer effective (as discussed above), and the DM density would be too 
high. The region of small A^ is thus not allowed beyond a narrow strip M^± > Ms- 

Also, for Ms = 75 GeV, if wc relax the constraints, we find solutions at even higher 
values of Mi. For example, if we disregard the AT constraint, or allow a somewhat larger 
value of r^DM, we also find solutions at Mi = 700 GeV, but then at a somewhat higher 
value of or A^,. 

For Ms ^ Mw, the annihilation SS — >■ W^W~ sets in, with a rate, which near 
threshold is controlled by 



(ax.)oc^. 



(7.1) 
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Figure 8: Allowed regions in the M^± — plane, for DM mass Ms = 110 GeV, with two values 
of Ma and lightest Higgs mass Mi = 200. 



with the final-state velocity 



As soon as this becomes sizable, the annihilation rate exceeds the value that is compatible 
with r^DM- Apart from the region discussed in Sec. 7.3, this annihilation mechanism thus 
provides an upper cut-off of the allowed region around Ms = 100 GeV. Eventually, as we 
will see in Sec. 8, for values of Ms <^ 550 GeV, the denominator in (7.1) will bring the 
annihilation rate down again, to an acceptable level. 

Below the cut-off, at Ms ^ 100 GeV, we again find an extension of the allowed range 
in m^, as Mi is increased towards 2Ms, where the annihilation via a neutral Higgs boson 
is resonant. This is illustrated in Fig. 7. For the particular case of Ms = 80 GeV, a wider 
range of solutions is found for Mi = 150 GeV than for Mi = 120 GeV. However, this has 
shrunk again at Mi = 155 GeV, and nothing is found for Mi = 160 GeV. In this region 
of Ms and Mi, the near-resonant annihilation via Hi provides too much depletion of the 
DM. 

7.3 Results for "new viable region" 

We confirm the existence of the new viable region identified by Lopez Honorez and Yamura 
for the IDM [52] . The existence of this region is due to cancellation between the four-point 
SSW~^W~ -coupling and the contribution via s-channel Hi exchange, proportional to X^. 
In our notation, the cancellation condition becomes [52] 



f^(M|-m2)=Ai=.--| 



or 



(7.3) 



m: 



= 3M| - iMf. (7.4) 

However, since we have adopted more tight constraints on OdM) we find a somewhat 
smaller allowed region. Examples are shown in Fig. 8, where we consider Ms = 110 GeV 
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and Ml = 200 GeV, with Ma = 205 GeV (left panel) and Ma = 215 GeV (right panel). 
The dominant loss mechanism here is SS — )• W~^W~ (four-point coupling, and rj^ t-channel 
exchange), together with a significant amount of SS — t- bb (via Hi s-channel exchange). 
We note a slight tendency for | | to increase with M^± , indicating that the s-channel 
increases in significance as the t-channel decreases. 

7.4 Summary of low-medium region 

For the low-to- medium mass region of Ms, an overview of allowed ranges of Mi is given in 
Fig. 9. The model is viable from about Ms ~ 6 GeV up to about 120 GeV. As a default 
value, we have taken Mi = 120 GeV, in the range favored by current electroweak SM 
fits [40]. The low range of Ms can be extended a bit by lowering Mi to 90 GeV. More 
interestingly, the value of Mi can be increased significantly, without any conflict with the 
data. It is typically restricted to Mi < 300 GeV, except for a narrow region around 
Ms ~ 75 GeV, for which values up to Mi ~ 600 GeV are acceptable. Beyond 75 GeV, 
it rapidly falls again (except for the region discussed in Sec. 7.3). This region around 
Ms = 75 GeV is a very interesting one, since the little hierarchy can be considerably 
alleviated [2], and, as we shall see in Sec. 10, the direct-detection cross sections are in this 
case within experimental reach. 

The cut-offs in Ms and Mi can be summarized as follows: 

• In the low range of Ms (~ 5 GeV), we can extend it a bit (lower Ms, higher Mi) by 
allowing a lower value of r^oM- 
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• In the high range of Ms 100 GeV), the cut-off is mainly due to too much DM 
annihilation in the Early Universe, via the gauge coupling. 

• The upper cut-off of Mi is mainly determined by the Odm and AT constraints. 
8. High DM Mass Regime 

We next study the model-parameter space when the DM is heavy. In analogy with the 
results for the simpler IDM [5], solutions are found for Ms ^ 545 GeV. In this mass 
range, annihilation via a single Higgs boson is not very efficient, whereas annihilations 
to two gauge bosons or two Higgs bosons are relevant. There is a lower cut-off around 
Ms ~ 540 — 550 GeV, below which the two-body annihilation to two gauge bosons, scaling 
like 1/M| (see Eq. (7.1)), is too fast to accommodate the observed value of 0dm- 

8.1 Scanning strategy 

In this high-mass region, the scanning is done differently from that of the low-medium-mass 
region. The main difference is that, for a viable model, we need Ma, M^± and m,, all to 
be close to Ms- There are actually three reasons for this: 

• A significant splitting would lead to "large" values of Aa, \ or Ac, and the Early- 
Universe annihilation to one or two Higgs bosons would be too fast, leaving too low 
a value for ^dm- 

• A significant splitting would require "large" values of Aa, A^ or Ac, and positivity or 
unitarity would be violated. 

• A significant splitting would lead to a value for AT in violation of the LEP data. 

The former constraints are stronger, such that the AT constraint has a negligible impact. 

Unless otherwise specified, in this section we consider Ma = Ms + 1 GeV. However, 
as we will discuss in Sec. 11, this splitting could be much smaller. 

8.2 Results 

Selected results on allowed regions in the M^±-m^ plane are displayed in Figs. 10-12. As 
indicated above, in this high-mass region, all inert-sector masses are rather degenerate, 
including the soft-mass parameter (implying that Ai is small). However, the way a 
correct DM density is obtained, is a bit different in the lower end of this high-mass region, 
where Ms ^ 550 GeV, from that of the higher-mass region, exemplified by Ms = 3000 GeV 
in Fig. 12. 

Figure 10 is devoted to Ms = 550 GeV. This is representative of the lower allowed 
value (in this high- mass region). Characteristic of this region is the annihilation via four- 
point gauge couplings, as illustrated by Eq. (4.5). These "loss" mechanisms involve not 
only SS annihilation, but also r]'^r]~ , AA, St]^ and Ai]"^ annihilations to gauge bosons. The 
cut-off around 540 — 550 GeV is due to the too high annihilation rate (at lower masses) [8] , 
scaling by Eq. (7.1). 
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Figure 10: Allowed regions (dark blue) in the Af^± — m,, plane, for DM mass Ms = 550 GeV, 
Ma = 551 GeV with lightest Higgs mass Mi = 120 GeV and 200 GeV. The thin sohd hne indicates 
m,j = Ms (and \l =0), whereas the dashed line gives M^± = Ma- 



As Ms is increased to higher values, the annihilation to gauge bosons will fall off, but 
annihilation to neutral and charged Higgs bosons start playing an important role, as shown 
by Eq. (4.6) for the case of Ms = 3000 GeV and Mi = 120 GeV. These annihilations are 
also due to four-point couplings, but originating from the potential, rather than being 
gauge couplings. Thus, their strengths can be adjusted by tuning the mass splittings in 
the inert sector. 

In Fig. 11 we display allowed regions in the M^-m^ plane, for Ms = 800 GeV. Com- 
pared to the case of Ms = 550 GeV, an allowed "ear"-shaped region has now developed. 
The interior is characterized by a high degree of degeneracy among the inert-sector mass 
parameters, which in turn leads to values of Al that are too small for an efficient annihila- 
tion of dark matter (or its partners) to final-state Higgs particles in the early universe, the 
value for would become too high. In fact, the dependence of J^dm on these couplings 
is quite steep: for the case {Ms, Mi) = (800, 120) GeV (left panel of Fig. 11) and the point 
{Mj^±,mrj) = {Ma, Ms), the value of J^dm is too high by a factor of two. On the other 
hand, outside the "ear", some A's will be large, the annihilation is too fast, and $7dm too 
small. 

As a final example, we show in Fig. 12 the cases of Ms = 3000 GeV and Ms = 
5000 GeV. Here, only the lower part of the "ear" is allowed. The upper part is forbidden, 
mostly due to the unitarity constraint. 
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Figure 11: Similar to Fig. 10, for Ms = 800 GcV. 



9. CP violation 



In order to illustrate the amount of CP violation which is available in the model we will 
plot so-called weak-basis-transformation invariants that have non-zero imaginary part. The 
benefit of studying invariants relies on the fact that they provide a measure of CP violation 
since any CP-violating observable must be a linear combination of the invariants (or their 
higher odd powers). The issue of the invariants was extensively discussed in the literature 
[54-59] in the context of the 2DHM. According to Gunion and Haber, there are three 
independent invariants Ji,2,3 that are sufficient to describe any CP-violating phenomenon 
which has its roots in the scalar potential (if Yukawa couplings are neglected). However, in 
the presence of three Higgs doublets the situation is much more involved and a complete 
study has not been performed up to date; for the existing attempts, see [54,55,59]. In 
order to discuss the invariants it is useful to write the potential in the following compact 
manner [56, 57] 

V{^i,^2,^3)=Y,-M^h + l:Z^lcMm^l^d), (9.1) 



where ^3 = r], and Y^i and Z^i^^ specify mass terms and quartic couplings, respectively. 
Formulae for Y^i and Z^i^^ in terms of the standard potential (2.3) parameters are given 
in Appendix B in (B.l). 

In general, for three doublets one should expect more invariants, however here for 
illustration, we discuss only Ji,2,3 defined in [58] for the two-doublet case. We generalize 
them in a very straightforward manner just by extending the range for indices from 1, 2 to 
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Figure 12: Similar to Fig. 10, for Ms = 3000 GeV and 5000 GeV, both with Mi = 120 GeV. 



1,2,3: 



h 
h 



J2 = vlvl^vlZ^eghZcfpfZeafd^aVdVhVr, 



(9.2) 



Of course, there exist additional invariants, however their determination is not necessary 
here. We have calculated the three invariants in the basis adopted in this paper (specified by 
the vacuum expectation values) for the general potential (2.3). The results are complicated 
quadratic and cubic polynomials in quartic coupling constants, shown in Appendix B, 
eqs. (B.3)-(B.5). However, if the dark democracy (2.9) is imposed the results simplify 
considerably: 



2„,2 



Im Ji 
Im J2 



i/(Ai-A2)ImA5 

A3 



(9.3) 



+ 



V" 



- ^4? - IA5H vf + 2(Ai - X2)Re\5vlvi 
((A2-A3-A4)'-|A5p)t;|] ImAs 
Ai)AcIm A5 



^ ^ ^(A3 + A4 + ReA 



^ ^ ^ (A3 + A4 + Re A5 - A2)AcIm A5 



(9.4) 
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Figure 13: Contour plots for averaged (over a's) imaginary parts of the invariants, 
Im Ji_2,3, illustrating the strength of CP violation, vs tan /J and M^, for [Ms, Ma, Mj^±,'mrj) = 
(75, 110, 90, 100) GeV, and (Mi, Ma, /i) = (120, 300, 200) GcV. 

2 2 

Im Js = ^(Ai - A2)(Ai + A2 + 2A4 + 2Ab)ImA5. (9.5) 

When finally we adopt the fact that Z2 is preserved by the vacuum (so = 0) then we 
obtain 

2 2 

ImJi = -^(Ai-A2)ImA5, (9.6) 
2 2 

Im J2 = [((Ai - A3 - A4)' - I Asp) vf + 2(Ai - A2)ReA5t;^2' 

- ((A2 - A3 - Xa? - IA5H vl] ImAs, (9.7) 

2 2 

Im J3 = ^(Ai - A2)(Ai + A2 + 2A4 + 2Ab)Im A5. (9.8) 

Note that, since Xa,b,c are real and the inert potential Vs{i]) is CP-conserving all the in- 
variants are proportional to the same CP-violating parameter Im A5 as in the case of the 
2HDM. In addition the dark democracy and (77) = imply that Im Ji^2 are identical to the 
corresponding invariants in the 2IIDM, whereas Im J3 differs by a term proportional to A;,, 
one of the quartic couplings between the non-inert and the inert sectors. 

For Ms = 75 GeV and -M^± = 90 GeV (and representative values for the other mass 
parameters), we show in Fig. 13 the imaginary parts of ^1,2,3 foi" parameters which are 
consistent with both experimental and theoretical constraints described in Sec. 3. Note 
that only low values of tan /3 are allowed. Here, averages over sets of a's are shown. It 
turns out that for the parameters adopted in Fig. 13, the contribution of A;, is small, of the 
order of ~ 5%. In the heavy DM case, where a high degree of degeneracy between inert 
masses is needed, the contribution from A^ is negligible. 

In Fig. 14 the corresponding predictions for the electron electric dipole moment are 
presented, both as an average over the same data set, and maximum (over a's) values, in 
units of [e 10~^^ cm], which is also the 1-a bound (see Eq. (3.10)). While typical values 
(left panel) are within 1 cr, the maxima (right panel) approach the cut-off, which is taken 
at 2 a. 
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Figure 14: Contour plots for averaged and maximal (over a's) electron electric dipole moment 
(in units [e 10^27 cm]), vs tan;9 and Af|, for {Ms,MA,M^±,m^) = (75,110,90,100) GeV, and 
(Mi,M2,/^) = (120,300,200) GeV. 



The presence of CP violation is also reflected in the distribution of 02 and 03, two of 
the three angles which determine the rotation matrix of the neutral Higgs sector, R. There 
are three limits of no CP violation, all identifiable in this plane [60]: 

Hi odd: 02 = =b7r/2, 03 arbitrary, 
H2 odd: 02 = 0, 03 = ±7r/2, 

H3 odd: a2 = 0, 03 = 0. (9.9) 

For the case M5 = 75 GeV, Mi = 120 GeV, we show in Fig. 15 the populated parts of 
this plane. There is a broad distribution of values, with no particular accumulation point. 
Thus, for the majority of these model points, CP is violated by a non-negligible amount 
(but still within the limits imposed by the EDM constraint). 



10. Direct detection 



The parameter points which give models compatible with particle-physics and DM con- 
straints, also give specific predictions for what signal should be observed in direct-detection 
experiments. In this regard, we have chosen to compare with the recent CDMS-II [61] and 
XENONIOO [62] results for spin-independent scattering. 

In Fig. 16 we compare the cross sections for our otherwise acceptable model points 
with these recent constraints. For each studied value of M5, we show a "column" of cross 
section values corresponding to different values of the other parameters. The cross section 
falls steeply with increasing mass M5. The different colors refer to the value of Mi, as 
specified in the caption. Values above the curves are excluded by the CDMS-II (solid curve) 
or the XENONIOO (dashed, red) experiments. We can here make a few observations: 
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a distributions 




Figure 15: Populated regions in the (02,0(3) plane, for Mg = 75 GeV, Mi = 120 GeV. 



• For a given value of Ms, the cross section tends to get lower for higher values of Mi 
(see the color coding). 

• At low masses, Ms ^ 10 GeV, the model is compatible with the bounds. 

• In the range 10 GeV ^ Ms ^ 50 GeV, most model points are excluded. 

• Around Ms ~ 60 — 80 GeV, much of the parameter space is again compatible with 
the bounds. 

• In the "high" region, the predicted cross sections are very low, implying that it would 
be difficult to test (exclude) the model in the near future. 

We should also stress that during the scanning, in the interest of covering as much of 
the "interesting" parameter space as possible, we did not evaluate the cross section. This 
is why some points violate the cross section constraint. 



11. LHC prospects 

At the LHC, one could imagine all the inert-sector scalars being pair-produced, 

pp SSX, AAX, SAX, Sr^-^X, Ar]^X, r)+r]-X, (11.1) 

followed by the decay of A or 77=*= to the lightest one, S. 

In favorable situations, decays involving -q^ could lead to observable signals. It was 
recently pointed out that in a related model [65], the combination of a small mass splitting 
between the charged scalar and the inert one, together with a small mixing angle, can lead 
to long-lived charged scalars that give displaced vertices in LHC detectors. In that model. 
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Cross sections 




Mg [GeV] 

Figure 16: Direct-detection cross sections compared with the CDMS-II (solid, [61]) and 
XENONIOO (dashed, [62]) bounds. Magenta: Mi < 120 GeV, green: 150 GeV < Mi < 230 GeV, 
blue: 300 GeV < Mi < 400 GeV, red: Mi > 500 GeV. 

the small mass splitting comes about from the assumption of unification at a high scale. In 
the model considered here, the small splitting is required by the appropriate prediction for 
the present DM abundance in the case of a heavy S (the DM candidate) . It is therefore of 
interest to check whether similar experimental signals are expected here as well. We split 
this discussion into two cases, according to the mass hierarchies. 

11.1 Ms < M^± < Ma 

The decay (via a virtual W) 

r/+ ^ Se+Ui (11.2) 

has several similarities to the familiar muon decay. The main differences are that (i) a 
scalar-scalar-vector vertex replaces a fermion-fermion-vector vertex, and (ii) one of the 
invisible final-state particles is massive. For the case of interest, M^± — Ms <^ M^±, the 
decay rate can be written as 

V = ^K±-^5)'. (11.3) 

We show in Fig. 17 mass differences, M^± — Ms, for a random subset of the allowed 
solutions found in Sees. 7 and 8. In addition, we show some more degenerate solutions 
(A^M = M^± — Ms = 0.1 GeV) for the high-mass region. In the low- medium-mass range, 
the mass splitting is obviously too large to give any interesting lifetime for rj^, as can also 
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Figure 17: Left: Mass differences, A*M = M^± — Ms vs Mg. Riglit: Corresponding 77* lifetime. 



be seen from Figs. 3-7. Even for M5 above the lowest allowed values of M^±, the co- 
annihilation Srj^ — )• W^'y prevents a mass splitting less than a few GeV. However, in the 
high-mass region, the mass splitting can be quite small, and longer lifetimes are possible. 

The experimental signature in this case would be the observation of an rj^ track from 
the production point up to the decay vertex (t]"*" — t- Si~^V£) followed by a kink corresponding 
to the decay and a track of the charged lepton. Although bosons of the inert doublet must 
be pair produced, such a kink does not depend on the nature of the accompanying boson 
being produced; be it r;^, yl or S, at least one kink is always there. Of course, in each case 
some missing energy (through the presence of SS in the final state) will also be helpful. 
The displacement of the decay vertex depends on the mass splitting and rj^ velocity, for 
example, for A^M = 0.1 GeV, the lifetime is 4.5 x 10~^ s, and with a velocity of /3 = 0.1, 
the decay length would be L = 13.6 m. This could easily be measured, and thus, e.g. for 
/? = 0.1 even a splitting as large as A^M = 0(1 GeV) could lead to observable effects. 

11.2 Ms < Ma< M^± 

In this case, if the mass splitting between A and S is small, — Ms <C M^, so any A 
produced would decay relatively slowly through A SZ* Sff with the width 



cos^ 0w 



{Ma - Ms)- 



(11.4) 



The experimental signature in this case would be a displacement between the production 
point and the A decay vertex, where the two tracks of the start. Again A would 
be accompanied by A, S or r/^, however at least one displaced vertex would be there, 
regardless what is the nature of the other inert boson. In each case missing energy will be 
carried away by two S bosons. In this case the displacement of the A decay vertex could 
be determined down to 0.2 cm, therefore the observation would be sensitive to values of 
Ma - Ms up to ~ 0.5 GeV. 
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12. Summciry 

We have explored in some detail the properties of the "IDM2" , an extension of the Inert 
Doublet Model ("IDM") [1,2] to the case of an additional non-inert doublet [7], allowing for 
CP violation. As compared with our earlier paper [7], we here implement the full positivity 
conditions, and run a full implementation of the model in micrOMEGAs. Allowed regions 
in the m,j-M^± and M5-M1 planes have been identified. Since our preliminary study [7] 
was restricted by < M^± , only the lower-right parts of Figs. 3-7 and Figs. 10-12 were 
accessible. We now see that, in most cases, a similar region above that diagonal (i.e., 
irijj > M^±) is also allowed. For some ranges of M5, it was found that the little hierarchy 
problem can be significantly reduced. 

In order to illustrate the strength of CP violation present in the model we generalize to 
the three-scalar-doublet case the three CP-sensitive weak-basis-transformation invariants 
•-^1,2,3 (defined originally for two scalar doublets). It turns out that Im Ji 2,3 ^ 0.5 — 3, 
which is five orders of magnitude more than the corresponding invariant in the Standard 
Model. 

Constraints on the model from direct-detection experiments have also been studied. 
In the region 10 GeV ^ Ms ^ 50 GeV, most model points are excluded by direct-detection 
experiments. 

Finally, if inert scalars are produced in vector-boson fusion at the LHC, we discuss 
possible lepton signatures from the decays of the heavier partners of the DM candidate. It 
turns out the measurable displacement of their decay vertex could provide an efficient way 
to test the model, at least for the heavy DM case. 
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Appendix A. Couplings of the inert sector 
Gauge couplings 

The quartic couplings of two gauge fields and two inert fields will have a "trivial" metric 
tensor factor, g^^'^, coupling the two gauge fields. In order to keep a "light" notation, we 
suppress this factor. Furthermore, wc denote the inert fields S, A and t]^. The heavy 
gauge fields are and Z, whereas the photon will be referred to as 7 (not A). In this 
notation, the quartic couplings involving two neutral inert fields are: 



SSW+W, AAW+W- : 




(A.la) 
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• 2 

SSZZ, AAZZ : ^ '\ ^ , (A. lb) 
whereas those involving two charged inert fields are: 

• 2 

rj+rj-W+W- : (A.2a) 

ri^rfll '■ 2ig^ SIT? 6w , (A. 2c) 

rf^'q'Z^: ig^ tauOw cos{2ew), (A.2d) 

(where for uniformity of notation we have not substituted e = ysin^^y) and finally the 
mixed ones are 

Sn^ZW^: -y^f^^ , (A.3a) 

' 2cos(9vF ' ^ ^ 

^.^^W^-: ^P4^, (A.3b) 
2 cos t^w^ 

. 2 

5r/±7T^T : ^sin^^y, (A.3c) 

Ar]^-iW^ : iysin^w'. (A.3d) 

There are no trilinear gauge couplings involving two identical inert fields. The non- 
vanishing trilinear gauge couplings (two inert fields and one gauge field) are: 

2^(^--^-)' (^-^^^ 

Srj^W^: ±^(ps-p±). (A.4b) 

Ar^^W^: |(pA-p±), (A.4c) 

r/V^: cos(2g^)(p+-p-), (A.4d) 
2 cos tfw 

V'^V~1 '■ — ig s^T^^wiP^ — P~), (A.4e) 

where all momenta ps, pA and p^ are incoming, and carry the Lorentz index of the vector 
field, in an obvious notation. 

Scalar couplings 

The scalar couplings can be expressed in a compact manner if we introduce the following 
notation: 

\l = \{K + \ + \c) = ^^^, (A.5a) 

\L=\{\a + \b-K)= "^^2 ' ( A-5b) 
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and 



Fj = cos PRji + sin PRj2, (A.6a) 

Fj = cos pRj2 - sin pUji . (A.6b) 

The latter quantities satisfy \Fj\ < 1, and \Fj\ < 1 since R is unitary. In particular, 
Fi = cos{(5 — ai) cos 02 . 

Trilinear couplings. The trilinear scalar couplings are: 



SSHj : 


- 2iXLvFj, 


(A.7a) 


AAHj : 


- 2i~XLvFj, 


(A.7b) 


SAHj : 


0, 


(A.7c) 


Sri'^H'^ : 


0, 


(A.7d) 


Ari^H"^ : 


0, 


(A.7e) 


V'^V'Hj : 


- iXaVFj, 


(A.7f) 


SSG^ : 


0, 


(A.7g) 


AAG^ : 


0, 


(A.7h) 


SAG° : 


- ivXc, 


(A.7i) 


Sri'^G^ : 




(A.7j) 


Ari^G^ : 


i \'"i^b - Ac). 


(A.7k) 



Prom Eqs. (A. 5) we note that the splittings M| — and — control the strengths of 
important trilinear couplings of two inert neutral fields to a Higgs field. Likewise, M^± — 
controls the strength of the inert charged fields to a Higgs field. If the dark democracy is 
lifted, all these couplings (A. 7) would be non-zero. 

Quadrilinear couplings. The quadrilinear ones involving two neutral inert fields are 



SSHjHj : - 2i{XL - XcR%), (A.8a) 

SSH^Hk-. 2iXcRj3Rk3, j^k, (A.8b) 

SSH+H- : - iXa, (A.8c) 

AAHjHj: -2iCXL + XcR%), (A.Sd) 

AAHjHk: -2iXcRj3Rk3, j k, (A.8e) 

AAH+H- : - iXa, (A.Sf) 

SAHjHj : - 2iXcRj3Fj, (A.8g) 

SAHjHk-. -iXc{FjRk3 + Rj3Fk), j ^ k, (A.Sh) 

SSG^G^ : - 2iXL, (A.8i) 

SSG+G- : - iXa, (A.Sj) 

SSH^G"^ : 0, (A.Sk) 
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SSG'^Hj : 0, (A.sr 

AAG^G'^ : -2iXL, {A.8m 

AAG+Q- : - iXa, (A.8n 

AAH'^G'^ : 0, (A.80; 

AAG^Hj : 0, (A.8p; 

SAG^G^ : 0, (A.Sq; 

SAG°Hj: -iXcFj, (A.Sr 

(A.8s; 

and those involving two charged inert fields are 

r]^r]^H^H^ : - 2iXc, {A.9a] 

rj+'q-HjHj : - iXa, {A.9h] 

r]+ri-HjHk: 0, j ^ k, (A.9c 

r]+r]-H+H-: - i(A„ + A^), (A.gd; 

r?±r/±G^G^ : - 2iXe, (A.de 

r]^r]^H^G^ : 0, {AM] 

n+'q-G'^G'^ : -iXa, (A.9g; 

r]+r]-G^Hj : 0, (A.9h; 

r]+r]-G+G- : -i{Xa + Xb), {A.9i 

ri+ri-H^G"^ : 0, (A.9j; 

and those involving one neutral and one charged charged inert field are 

Srj^H^Hj : ^{Xb + Xc)Fj ± ^(A^ - Xc)Rj3, (A.lOa; 

Ari^H^Hj: ^(X^ ~ X,)Fj - '-{Xb + Xc)Rj3, (A.lOb; 

5?7±G^G0 : ± ^{Xb - Ae), (A.lOc; 

Sr]^H'^G^ : 0, (A.lOd; 

Srj^G^Hj : - ^{Xb + X,)Fj, {AAOe 

Arj^G^G'': -'-(Xb + Xe), (A.lOf 

Ari^H'^G'^ : 0, (A.lOg; 

ArfG^Hj-. T\{\h-K)Fy (A.lOh; 



In Eq. (A.7)-(A.10), the couplings listed as zero, would be non-zero if the dark democracy 
is lifted. On the other hand, couplings not listed, are absent also when the dark democracy 
is lifted. 
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Couplings involving four fields, all from the inert doublet, are 



SSSS : 


— SiXfj, 


AAAA : 


— SiXjj, 


SSAA: 


- iXr,, 


SSri+rj- : 


- iXr,, 


AAri'^rj~ : 


- iXr,, 




— 2zA^. 



(A.lla) 
(A.llb) 
(A.llc) 
(A.lld) 
(A.lle) 
(A.llf) 



Appendix B. Invariants Im Ji,2,3 



The tensors Z^^^d adopted in (9.1) can be expressed through the standard parameters used 
in (2.3) as follows: 





■^llll 


= Ai, 




■^2222 


= A2, 




■^3333 


= A^, 


-^1122 


= -^2211 


= A3, 


-^1221 


= -^2112 


= A4, 


-^1133 


= -^3311 


= A1133 


■^2233 


= -^3322 


= A2233) 


■^1331 


= -^3113 


= A1331, 


■^2332 


= -^3223 


= A2332 




■2^1212 


= As, 




■^2121 


— \* 




-^1313 


= A1313 




■^3131 


= A^3i3, 




■^2323 


= A2323) 




■^3232 


— \* 

— ^^2323 



(B.l) 
We also have 

vi = vi/v, V2 = V2/V, V3 = V3/V. (B.2) 
The invariants Im Ji,2,3 expressed in the basis adopted here then read: 

Im Ji = -;^(Ai - A2 + A1331 - A2332)('yilmA5 + i;|lmAi3i3), (B.3) 



Im J2 = ^ (f |lm A5 + v^ImX 



1313 J 



((Ai - A3 - A4)' - lAsP) vf + 2(Ai - A2)ReA5i;^i - ((A2 - A3 - A4)' - IA5P) vl 

+ 2t'l^3 [(A1133 + A1331 — A2233 — A2332)(Aii33 + A1331 — A3 — A4 + Re A1313) 

-ReAi3i3(A3 + A4 - Ai) - ReA2323ReA5] 

+ 2^2 1"! [(A1133 + A1331 — A2233 — A2332)(A2233 + A2332 — A3 — A4 + Re A2323) 

+Re A2323(A3 + A4 - A2) + Re Ai3i3Re A5 + Im Ai3i3lm A5] 

+ [|Ai3i3p - (A1133 + Ai33i)^ 

+ (A2233 + A2332 — Re A2323)(A2233 + A2332 + ReA2323) 
+2(Aii33 + A1331 — A2233 — A2332)A,j] 

(B.4) 



2vfvL . . . vU^ 



^ 3lmAi3i3lmA5 - ^(ImAi3i3)^ 



,2 



Im J3 = ^(Ai - A2 + A1331 - A2332)(Ai + A2 + 2A4 + A1331 + A2332)(t'ilm A5 + i;|lm A1313). 

(B.5) 
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where, in order to retain generality we kept f3/\/2 = (77) 7^ 0, and made use of the relation 

vflm A1313 + i;|lm A2323 = 0, (B.6) 



which emerges from the minimization conditions for ^3 7^ 0. 

In the dark democracy case, these results simplify to those given in Eqs. (9.3)-(9.5). 
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